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Abstract — The  identification  and  discrimination  of  unexploded 
ordnance  using  low-frequency  electromagnetic  induction  is  an 
expensive  and  difficult  process,  typically  beset  by  low  data  diver¬ 
sity  and  high  positioning  uncertainty.  In  this  paper,  we  present 
the  Man-Portable  Vector  (MPV)  sensor,  a  new  time-domain  in¬ 
strument  designed  to  remedy  these  shortcomings  by  measuring 
all  three  vector  components  of  the  secondary  magnetic  field  at 
five  distinct  points  around  each  transmitter  location.  The  MPV 
also  has  a  laser  positioning  system  that  can  give  its  location  with 
millimeter  precision.  After  describing  the  instrument  in  detail,  we 
study  its  performance  in  various  sets  of  measurements,  using  the 
tensor  dipole  model  to  analyze  the  data.  We  find  that  the  sensor 
can  detect  deeply  buried  targets  and  identify  some  standard  ord¬ 
nance  items.  It  can  also  resolve  separate  targets  in  cases  where  two 
objects  share  the  field  of  view  and  produce  overlapping  signals.  A 
new  incarnation  of  the  MPV,  the  MPV-II,  is  in  an  advanced  stage 
of  development. 

Index  Terms — Electromagnetic  induction  (EMI),  man-portable 
vector  (MPV)  sensor,  unexploded  ordnance  (UXO). 


I.  Introduction 

THE  PRESENCE  of  buried  unexploded  ordnance  (UXO) 
in  former  battlefields  and  decommissioned  firing  ranges 
is  an  increasingly  serious  humanitarian  and  environmental 
problem  worldwide  [1],  [2],  with  timescales  measured  in 
decades  and  yearly  casualty  rates  in  the  hundreds.  UXO 
cleanup  is  extremely  expensive  and  laborious,  reaching  some 
$1.4  million  per  acre  in  the  United  States  [3],  because  most 
environments  that  contain  UXO  are  also  cluttered  with  harm¬ 
less  items  that  get  detected  just  as  readily  and,  in  the  absence 
of  further  information,  must  be  treated  as  dangerous.  Thus  the 
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Fig.  L  MPV  sensor  with  the  ArcSecond  positioning  system.  The  sensor  head 
is  about  75  cm  in  diameter;  see  Fig.  2  for  the  other  dimensions. 

technical  problem  is  not  one  of  detection  but  one  of  identifi¬ 
cation  and  discrimination.  Electromagnetic  induction  (EMI) 
sensing,  in  which  a  primary  field  established  by  a  sensor  in¬ 
duces  eddy  currents  and  magnetic  response  in  metallic  bodies 
[4]  in  a  frequency  range  at  which  the  ground  is  essentially 
transparent,  is  considered  a  most  promising  technique  in  the 
field;  however,  it  is  not  exempt  from  difficulties. 

Induction  sensing  occurs  at  very  low  electromagnetic  fre¬ 
quencies,  which  means  that  one  must  work  at  very  low  reso¬ 
lutions  and  forfeit  the  possibility  of  imaging  to  search  instead 
for  patterns  or  signatures  in  the  data.  This  shortage  of  infor¬ 
mation  is  exacerbated  by  the  very  small  spatial  diversity  of  the 
measurements,  which  often  consist  of  a  single  field  component 
measured  at  a  single  altitude  at  only  a  few  points,  and  by  the 
strong  influence  of  sensor  positioning  and  its  uncertainty  on  the 
detected  signal.  New-generation  instruments  like  the  Berkeley 
UXO  Discriminator  (BUD)  [5],  [6],  TEMTADS  [7],  and  the 
MetalMapper  [8],  [9]  attempt  to  ameliorate  this  situation  by  pro¬ 
viding  arrays  of  transmitters  and  receivers  in  precisely  known 
configurations;  these  systems,  however,  are  large  and  heavy  and 
cannot  be  readily  used  in  rough  or  treed  terrain  that  does  not 
allow  vehicular  access. 

The  Man-Portable  Vector  (MPV)  sensor  [10],  shown  in  Fig.  1 
and  schematically  in  Fig.  2,  has  been  designed  in  response  to  the 
previous  concerns.  Developed  by  G&G  Sciences,  the  MPV  is 
a  mono/multistatic  ultrawideband  time-domain  handheld  EMI 
instrument  that  measures  all  three  components  of  the  secondary 
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Fig.  2.  Sketch  of  the  MPV  from  behind  (top  panel)  and  from  above  (bottom 
panel).  The  red  circles  of  radius  37.5  cm  are  the  two  transmitter  coils,  described 
more  fully  in  Section  II- A.  Their  finite  heights  have  been  neglected.  The  black 
squares  of  side  10  cm  represent  the  receivers,  with  R1  at  the  center,  R0  located 
30.6  cm  above  it  (in  the  +zMpv  direction),  and  the  others  separated  from  R1  by 
39.3  cm:  R2  in  the  — xMPV  direction,  R3  in  +yMPV,  and  R4  in  +xMPV.  These 
are  described  in  detail  in  Section  II.B.  The  ArcSecond  receivers  of  Section  II-C, 
arranged  in  an  isosceles  right  triangle  above  and  behind  the  receiver  assembly, 
are  represented  by  large  blue  dots. 


field  at  five  different  locations — in  particular  at  two  different 
heights — with  user-adjustable  temporal  resolution  and  data-ac- 
quisition  window,  either  in  static  or  dynamic  measurement  con¬ 
figurations.  It  is  equipped  with  a  laser  positioning  system  that 
tracks  its  location  with  sub-centimeter  precision.  It  also  records 
the  transmitter  current  at  every  point,  which  can  be  used  to 
cancel  errors  in  this  current  due  to  instrument  drift.  These  fea¬ 
tures,  along  with  its  portability,  make  the  MPV  a  highly  versatile 
instrument  that  provides  a  wealth  of  high-quality  information. 

Several  previous  studies  [1 1]— [1 3]  have  shown  that  the  MPV 
is  indeed  a  capable  tool  for  detection  and  identification  of  buried 
items.  MPV-collected  data  have  been  analyzed  using  the  dipole 
model  [14],  the  Standardized  Excitation  Approach  [15]  and  a 
generalized  version  thereof  [12],  and  the  Normalized  Surface 
Magnetic  Source  (NSMS)  method  [11],  [16],  [17]  to  detect  and 
identify  single-  and  multiple-target  UXO  scenarios  in  the  labo¬ 
ratory  and  at  various  UXO  sites.  Future  plans  include  the  adap¬ 
tation  to  the  MPV  of  the  Orthogonal  Volume  Magnetic  Source 
model  [18]  and  of  spheroidal  [19],  [20]  and  ellipsoidal  [21]  an¬ 
alytic  procedures. 

In  this  paper,  we  aim  to  present  the  MPV  in  detail  and  use 
MPV  data  to  perform  some  examples  of  identification  tasks 
using  the  tensor  dipole  model,  which  has  worked  well  with  other 
instruments,  both  in  the  frequency  domain  [22]-[26]  and  in  the 
time  domain  [27].  This  method,  which  characterizes  a  target 
using  a  symmetric  polarizability  tensor  [28],  allows  reasonable 


estimation  of  the  location  and  orientation  of  a  buried  object  and 
yields  an  intrinsic  electromagnetic  signature  that  can  be  ana¬ 
lyzed  further — for  example,  by  checking  it  against  a  standard¬ 
ized  set  previously  stored  in  a  library.  The  next  step — use  of  this 
identifying  information  to  classify  buried  targets  and  ultimately 
to  discriminate  UXO  from  clutter — is  beyond  the  scope  of  this 
paper  but  has  been  performed  successfully  starting  from  MPV 
data  [15],  [16]. 

The  paper  is  organized  as  follows.  In  Section  II,  we  present 
the  MPV  sensor,  describing  its  features  and  positioning  system 
and  addressing  some  of  the  special  data-processing  procedures 
required  by  them.  In  Section  III,  we  introduce  the  dipole 
model  and  use  it  to  give  a  few  examples  of  one-  and  two-target 
identification  and  discrimination  starting  from  MPV  data.  In 
Section  IV,  we  present  the  MPV-II,  a  forthcoming  streamlined 
and  more  ergonomic  version  of  the  MPV,  and  in  Section  V  we 
conclude. 


II.  MPV  Sensor 

The  MPV  is  a  hybrid  monostatic/multistatic  sensor.  The 
transmitters  and  receivers  are  rigidly  attached  and  move  in 
lockstep,  but  for  each  sensor  location  (transmitter  sequence 
or  “data  shot”)  there  are  five  receiver  locations  and  three 
vector  components  at  each  receiver,  for  a  total  of  15  usable 
measurements.  Sensor  positioning  is  currently  tracked  using 
an  ArcSecond  laser  system  (the  gray  PVC  masts  appearing 
in  Fig.  1)  that  maintains  accuracies  on  the  order  of  a  few 
millimeters.  The  sensor  is  designed  to  achieve  a  high  SNR  with 
a  large  dynamic  range.  It  can  in  principle  be  used  for  ordinary 
surveying  or  for  close  interrogation  of  anomalies  previously 
detected  by  other  methods. 

The  MPV  can  be  operated  in  two  distinct  data-acquisition 
modes:  static  and  dynamic.  In  both  modes,  the  MPV  collects  po¬ 
sitioning  and  magnetic  field  data  at  a  rate  of  10  Hz.  The  100-ms 
cycle  is  described  further  in  Section  II-A.  In  static  data-acquisi- 
tion  mode  the  instrument  is  placed  in  a  fixed  location — typically 
at  a  set  of  approximate  grid  points  over  a  specific  anomaly — and 
the  10-Hz  data  are  stacked,  or  averaged,  for  an  amount  of  time 
that  can  be  selected  by  the  user.  This  mode  results  in  a  higher 
SNR  compared  to  that  of  the  dynamic  data  mode  due  to  the 
noise-reducing  effect  of  integration,  but  the  resulting  data  are 
fairly  sparse.  Using  the  MPV  in  dynamic  mode  results  in  10-Hz 
data  collection  points  as  the  MPV  is  freely  waved  around  above 
an  anomaly.  While  the  SNR  is  lower  in  this  mode,  the  data 
are  potentially  more  diverse  because  the  instrument  can  be  ma¬ 
neuvered  to  interrogate  the  anomaly  from  different  angles  and 
depths.  The  data  are  also  more  abundant:  half  a  minute’s  worth 
of  target  illumination  results  in  300  x  5  x  3  =  4500  data  points, 
not  including  the  number  of  time  gates  chosen  for  each.  We 
shall  concentrate  on  the  static  mode  in  this  paper;  we  note,  how¬ 
ever,  that  it  is  equally  possible  to  perform  UXO  discrimination 
starting  from  dynamic  data  [11],  [29]. 

A.  Transmitter  Coils 

The  transmitter  assembly  of  the  MPV  consists  of  two  circular 
parallel  coils  of  radius  37.5  cm  whose  centers  lie  on  a  line  per¬ 
pendicular  to  the  planes  of  the  circles;  the  distance  between  the 
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Fig.  3.  MPV  transmitter  cycle.  The  plot  shows  a  typical  bipolar  on/off  cycle 
of  the  MPV  transmitter  assembly.  The  total  duration  is  100  ms,  divided  into 
four  stages  of  equal  duration  (separated  by  vertical  dotted  lines):  A  step-on 
that  plateaus  at  about  4.5  A,  a  step-off  (shown  in  detail  in  the  top  left  inset),  a 
step-on  of  opposite  polarity,  and  another  step-off.  The  bottom  right  inset  shows 
the  second  (“negative”)  turn-on  stage;  we  use  a  logarithmic  y-axis  to  bring  at¬ 
tention  to  a  small  peak  that  occurs  between  the  first,  more  abrupt  rise  and  the 
damped  one  that  ends  in  a  stable  value.  In  actual  measurements  only  the  step-on 
currents  are  recorded.  The  data  acquisition  system  has  a  sampling  rate  of  2  ^s; 
the  plot  thus  contains  5  x  104  data  points.  Each  point  has  resulted  from  stacking 
100  measurements. 

bottom  of  the  upper  transmitter  and  the  top  of  the  lower  trans¬ 
mitter  is  1 1 .6  cm.  The  bottom  of  lower  transmitter  coil  is  0.7  cm 
from  the  bottom  surface  of  the  plywood  framework  that  houses 
it.  Each  coil  is  about  4.8  cm  tall  and  is  composed  of  14  helically 
wound  copper- wire  loops,  each  of  which  has  diameter  2  mm 
and  is  vertically  separated  from  its  neighbors  by  approximately 
1.5  mm  end  to  end.  These  coils  vertically  bracket  most  of  the 
receiver  assembly  to  be  discussed  in  the  next  subsection. 

The  transmitter  works  in  the  bipolar  on/off  cycle  shown  in 
Fig.  3.  (These  settings  are  in  fact  adjustable;  here  we  restrict  our 
attention  to  the  numbers  used  most  often  in  the  measurements 
reported  below.)  A  gradual  but  fairly  rapid  buildup  of  current 
occurs  at  the  beginning  of  each  cycle;  after  a  couple  of  millisec¬ 
onds  the  current  plateaus  and  stabilizes  at  about  4.5  A.  At  25  ms 
it  is  shut  off  abruptly  and  kept  at  zero  for  25  ms  more.  Midway 
through  the  cycle  the  current  is  built  up  again,  with  the  polarity 
reversed  in  order  to  limit  the  magnetization  of  the  target  during 
data  collection  and  suppress  long-lived  eddy  currents,  and  then 
shut  off  again  at  75  ms  for  another  25  ms  that  conclude  the  cycle. 
Sharp  spikes  follow  the  shutoffs  and,  to  a  lesser  extent,  the  onsets 
of  the  buildup  stages;  these  are  highlighted  in  the  figure  insets. 

To  make  the  data  files  more  compact,  only  the  current  through 
the  transmitter  coils  is  measured  during  transmitter-on  time; 
when  the  transmitter  is  turned  off,  only  the  voltages  across  the 
15  receiver  coils  are  recorded.  The  MPV  thus  gives  a  picture 
of  the  transmitter  current  when  it  is  on,  at  different  times  than 
the  receiver  signals  but  sampled  identically.  The  reason  for  sam¬ 
pling  the  transmitter  current  is  to  measure  it  near  the  end  of  its 
cycle,  when  it  reaches  its  peak,  because  that  is  the  magnitude 
of  the  driving  signal  that  energizes  a  target:  the  driving  field  at 
the  target  ( dB / dt  from  the  transmitter  current)  is  approximately 
an  impulse  with  a  magnitude  proportional  to  the  current  in  the 
transmitting  loop  just  before  it  turns  off.  The  data  have  to  be 


Fig.  4.  Close-up  view  of  one  of  the  MPV  cubes  (R0  in  this  case)  showing  the 
three  receiver  coils.  The  numbers  in  the  tape  measure  increase  in  the  -%pv 
direction. 


normalized  with  respect  to  the  transmitter  current  because  the 
latter  varies  noticeably  from  point  to  point.  Usually  the  mea¬ 
sured  voltage  at  each  receiver  coil  is  divided  by  the  transmitted 
current  read  at  the  last  time  sample. 

In  the  modeling  examples  that  follow  we  neglect  the  heights 
of  the  coil  assemblies  and  simulate  the  transmitters  using  an  an¬ 
alytic  expression  for  the  primary  field  due  to  an  infinitesimally 
thin  circular  loop  of  radius  x  placed  flat  on  the  x-y  plane  and 
centered  at  the  origin  [30],  [31]: 


Hpr(r)  =  -  V  x 

7 r 


K(k)  -  E(k) 


0) 


in  cylindrical  coordinates,  where  the  modulus  of  the  complete 
elliptic  integrals  K  and  E  [32],  [33]  is 


k2 


4  PX 

(X  +  P )2  +  ’ 


(2) 


We  take  the  two  thin  loops  to  be  separated  by  15.6  cm. 


B.  Receiver  Cubes 

The  MPV  receiver  assembly  consists  of  five  triaxial  coils 
wound  around  plastic  cubes,  numbered  R0  through  R4.  Fig.  2 
shows  a  diagram  of  the  sensor  as  seen  from  above  and  from  be¬ 
hind.  In  what  follows  we  refer  to  the  coordinate  system  shown 
in  the  figure  as  the  “MPV”:  the  origin  is  at  the  center  of  receiver 
Rl;  the  rr-axis  points  from  the  origin  to  the  center  of  receiver 
R4,  the  y-axis  points  to  the  center  of  receiver  R3,  and  the  z-axis 
points  to  the  center  of  receiver  R0.  In  that  system,  cube  R0  is 
located  above  center  (z  =  30.6  cm),  cube  R2  is  left  of  center 
(x  =  —39.3  cm),  cube  R3  is  forward  of  center  ( y  =  39.3  cm), 
and  cube  R4  is  right  of  center  ( x  =  39.3  cm). 

Fig.  4  shows  a  close-up  view  of  one  of  the  cubes.  Each 
air-cored  receiver  coil  is  wound  helically  on  a  set  of  ten 
grooves  etched  on  the  surface  of  the  cube;  36-gauge  wire  is 
used  throughout.  Each  groove  is  1 .2  mm  wide  and  is  separated 
from  its  neighbors  by  2.0  mm;  thus,  each  set  of  grooves  has  a 
total  width  of  0.12  x  10  +  0.2  x  9  =  3.0  cm.  From  Faraday’s 
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law,  the  measured  signal  must  be  interpreted  as  the  surface 
integral  of  the  time  derivative  of  the  secondary  magnetic  flux 
density;  the  dimensions  of  the  MPV  receiver  coils  have  been 
chosen  so  the  signals  are  as  similar  loop-to-loop  as  possible. 
The  grooves  corresponding  to  the  z-coil  (the  horizontal  one 
in  Fig.  4)  are  1.0  mm  deep;  each  groove  has  12  turns  wound 
around  it,  giving  the  coil  an  effective  flux-measuring  area 
of  120  x  (10.0  -  0.1  x  2 )2  =  11524.8  cm2.  The  a;-coil, 
perpendicular  to  the  tape  measure  in  Fig.  4,  is  wound  around 
4.8-mm  deep  horizontal  grooves  and  1.0-mm  deep  vertical 
ones  and  contains  13  turns  per  groove;  its  effective  area  is 
130  x  (10.0  -  2  x  0.48)  x  (10.0  -  2  x  0.1)  =  11517.0  cm2. 
Finally,  the  y- coil  surrounds  4.66-mm  deep  grooves,  and 
has  14  turns  per  groove,  which  gives  it  an  effective  area  of 
140  x  (10.0  —  2  x  0.466)2  =  11512.0  cm2.  In  the  inversions 
of  Section  III  we  take  all  receivers  to  be  10  x  10  cm2;  we 
also  implicitly  divide  out  the  number  of  loops  and  neglect  the 
groove-to-groove  separation.  These  effects  are  not  too  dramatic 
unless  the  target  is  very  close  to  the  sensor,  as  has  also  been 
found  to  happen  with  other  EMI  instruments  [34] . 

The  MPV  records  a  user- selectable  number  of  time  gates 
ranging  from  about  60  [is  to  typically  25  ms.  Common  settings 
include  30  or  35  logarithmically  spaced  gates,  but  the  MPV  can 
also  sample  more  than  100  time  channels  (as  in  Section  III-B1). 
The  signals  from  the  cubes  are  amplified  and  sampled  in  a  data- 
acquisition  system  tethered  to  the  sensor  (the  backpack  in  Fig.  1). 


C.  ArcSecond  Positioning 

Access  to  the  full  3D  response  of  a  target  makes  it  necessary 
to  have  a  full  3D  positioning  system  that  accounts  for  the  loca¬ 
tion  and  tilt  angles  of  the  sensor  at  any  point.  Some  existing  EMI 
sensors  have  used  GPS  for  this  purpose  [5],  [6],  [35]  though  the 
available  resolution,  on  the  order  of  centimeters,  serves  to  record 
anomaly  locations  for  further  study  but  sometimes  does  not  suf¬ 
fice  to  support  discrimination. 

The  MPV  is  equipped  with  an  ArcSecond  positioning  system 
that  consists  of  three  receivers,  SI,  S2,  and  S3,  whose  locations 
are  known  with  respect  to  a  coordinate  system  (the  “ASI”)  pro¬ 
vided  by  the  positioning  device.  The  ArcSecond  receivers  are 
arranged  in  an  isosceles  right  triangle  of  side  25  cm  placed  so 
that  the  line  pointing  from  the  centroid  of  the  triangle  to  the  right 
angle  coincides  with  yMpv  and  the  hypotenuse  coincides  with 
xmpv-  The  centroid  of  the  triangle  is  located  22.4  cm  behind 
receiver  R0  (see  Fig.  2.)  In  terms  of  the  ASI  receiver  locations 


we  have  the  MPV  unit  vectors  [36] 

rS2  -  rS3 

^MPV  =  ,  , 

|rS2  -  rS3| 

(3a) 

_  rsi  —  |(rS2  +  rS3) 

|rsi  -  KrS2 +rS3)| 

(3b) 

zmpv  =  ^mpv  x  i/MPV 

(3c) 

which  in  turn  yield  the  rotation  matrix  (also  called  “direction 
cosine  matrix,  or  DCM”) 


R  = 


xmpv  •  x 
2/mpv  *  -T’ 
zmpv  •  x 


^mpv  •  y  ^mpv  •  z 
Vmpv  *  V  2/mpv  *  ^ 
zmpv  •  y  %pv  *  z 


(4) 


At  each  step  the  DCM  can  be  polished  (i.e.,  made  “more  orthog¬ 
onal”)  by  iterating  [37] 

Rn+1  =  i(Rn  +  (RT)“1)  (5) 

starting  from  the  initially  measured  Rq.  Usually  three  iterations 
are  sufficient.  From  the  DCM  we  can  compute  the  Euler  angles 


in  the  yaw-pitch-roll  convention  [38],  modified  so  that  x 

->■  ?/• 

y  — ►  x,  z  — »  —z  [36]: 

/  J  2/mpv  *  & 

(p  =  arctan  - - — - 

Yaw 

(6a) 

2/mpv  •  y 

6  =  arcsin  ?/MPV  •  z 

Pitch 

(6b) 

,  ,  xmpv  •  z 

—%)  =  arctan  - - - 

Roll. 

(6c) 

^MPV  ‘  Z 

A  relevant  measurable  parameter  is  the  “tilt  angle”  of  the  MPV, 
defined  by 


7  =  arccos  (zMpv  *  z)  (7) 

which  in  terms  of  the  Euler  angles  is 

cos  7  =  cos  cos  6.  (8) 

This  expression  is  independent  of  yaw  and  has  a  different  inter¬ 
pretation  for  cases  with  pure  pitch  than  in  cases  with  pure  roll, 
allowing  us  to  check  the  correctness  of  each  angle  separately. 

Next  we  locate  the  MPV  and  its  receivers.  The  ASI  triangle 
is  centered  at  the  point 

Rc  =  ^(rsi  +  rS2  +  rss)  (9) 

and  from  the  information  in  Fig.  2  we  can  locate  every  other 
point  on  the  sensor.  For  example,  the  “center”  of  the  MPV  (i.e., 
the  position  of  receiver  Rl)  is 


rRi  =  Rc  +  22.4 t/MPV  —  ((81.5  —  5.9/2)  +  21.8/2)  zmpv cm. 

(10) 

Fig.  5  shows  an  example  measurement.  With  the  MPV  ini¬ 
tially  flat,  the  experimenter  lifted  the  right-hand  side  (xmpv)  of 
the  sensor  to  an  angle  of  about  20  degrees  and  then,  always  with 
a  point  on  the  bottom  rim  of  the  instrument  in  contact  with  a  flat 
surface,  and  the  opposite  side  in  the  air  at  about  the  same  initial 
angle,  in  one  smooth  motion  performed  the  cycle 

xmpv  up  — >  2/Mpv  up  — >  -xMpv  up 
*  — 2/mpv  UP  *  *^mpv  up 

and  at  the  end  laid  the  MPV  flat  again,  pointing  in  the  initial  di¬ 
rection.  (This  motion  somewhat  resembles  that  of  a  slowly  spin¬ 
ning  coin.)  Fig.  5(a)  shows  the  location  of  the  ASI  triangle  as 
the  sensor  moved,  taken  directly  from  the  MPV  readings.  The 
protruding  lines  (not  to  scale)  represent  the  unit  vector  zmpv, 
which,  as  expected,  is  tilted  at  a  constant  angle  and  is  seen  to 
trace  a  circle.  The  dots  at  the  bottom  show  the  receiver  Rl, 
which  barely  moves  during  the  process.  Fig.  5(b)  depicts  the 
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Fig.  5.  Test  of  the  ASI  positioning  system  on  the  MPV.  An  experiment  was  carried  out  in  which  the  MPV  was  tilted  and  cycled  in  a  smooth  motion  reminiscent 
of  that  of  a  tilted  spinning  coin,  (a)  Location  of  the  ASI  triangle  (brown)  as  the  sensor  moves,  the  footprint  of  the  (barely  moving)  bottom  central  receiver  R1  (red), 
and  the  tilted  zMPV-axis  (blue  lines,  with  arrows  removed  to  decrease  clutter).  The  sizes  of  the  triangle  and  a  receiver  are  shown  for  reference,  (b)  Euler  angles 
(dashed  lines)  in  the  yaw-pitch-roll  convention,  showing  good  agreement  between  the  two  expressions  (7)  (solid  line)  and  (10)  (circles)  for  the  tilt  angle  7. 


Euler  angles  as  computed  from  (6);  the  agreement  between  ex¬ 
pressions  (7)  and  (8)  for  the  tilt  angle  7  shows  that  the  pitch, 
yaw,  and  roll  are  determined  correctly. 

III.  Example  Inversions  Using  the  MPV 

A.  Dipole  Model 

A  popular  and  useful  technique  for  UXO  discrimination  with 
EMI  sensors  models  a  buried  target  by  a  point  dipole.  The  model 
has  been  implemented  in  the  frequency  domain  [24],  [22]  and 
in  the  time  domain  [27]  and  has  been  found  to  fit  measured  data 
adequately  when  the  targets  are  small  enough — or  far  enough 
away  from  the  sensor — that  the  primary  field  can  be  assumed 
uniform  along  their  extent,  and  when  their  composition  is  ho¬ 
mogeneous  enough  that  there  are  no  conflicting  signals  from 
different  parts  of  the  object  [23],  [39]. 

The  dipole  model  postulates  that  for  a  given  object  it  is  pos¬ 
sible  to  find  a  set  of  orthonormal  axes  (a,6,c)  such  that  an 
external  monochromatic  magnetic  field  pointing  along  one  of 
those  directions  induces  a  steady-state  dipole  moment  parallel 
to  it: 

ma  =  paHV\  mb  =  pbHl\  mc  =  pcH?.  (11) 

For  a  primary  field  pointing  in  an  arbitrary  direction  in  the  lab¬ 
oratory  system,  (11)  becomes  [23] 

m  =  Ratc  A  Rj6c  Hpr  =  B  Hpr  (12) 

where  Ra^c  is  an  orthogonal  DCM  like  (4)  and  A  is  a  diag¬ 
onal  matrix  whose  elements  are  the  from  (11).  The  “polar¬ 
izability  tensor”  B  is  a  symmetric  3x3  matrix  with  dimensions 
of  volume  [28].  The  secondary  field  is  given  by  the  standard 
dipole  expression  [40],  [41] 

1 

4nR 5 


Fig.  6.  Grid  used  for  MPV  test-stand  measurements.  The  data  runs  start  at  point 
ri  on  the  upper  left  corner  and  follow  the  blue  line  until  point  r89  =  190x  + 
80y  cm  on  the  lower  right  corner;  each  is  preceded  and  followed  by  calibration 
measurements,  some  1 80  cm  away  from  the  center.  Superimposed  is  a  schematic 
diagram  of  the  sensor,  shown  in  more  detail  in  Fig.  2. 

where  R  =  r  —  r',  r  is  the  observation  point,  and  r'  is  the 
location  of  the  point  dipole.  In  observance  of  Faraday’s  law, 
EMI  receivers  measure  not  the  field  (13)  but  rather  (the  negative 
of)  the  time  derivative  of  the  flux  through  them.  The  time-do- 
main  signal  measured  at  the  kth  time  gate  by  receiver  Rp  (p  = 
0, . . . ,  4)  along  the  uth  (y  =  z,y,x)  direction  when  the  sensor 
is  centered  at  17  is  thus 

&  nv  (j^j  5 

jRn(rj) 

■  -A_(3R(r  •  p(r j,tk))  -  R2p(rj,tk)) 


Hsc(r)  = 


(3  R(R  •  m)  —  R2  m) 


(13) 


(14) 
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Fig.  7.  Raw  data  taken  by  the  MPV  in  Vicksburg.  This  particular  data  set  has  an  upright  81-mm  projectile  as  target,  located  45  cm  below  the  sensor,  (a)  Point-to- 
point  values  at  the  10th  time  gate  (0.168  ms  after  shutoff  in  this  case)  in  the  (from  top  to  bottom)  z,  x,  and  y  directions,  (b)  Time  development  (in  absolute  value) 
at  point  r45,  right  at  the  center  of  the  grid.  The  data  has  been  normalized  by  the  current  at  every  point,  and  a  background  field  taken  some  180  cm  away  with  no 
targets  present  (and  smaller  by  four  orders  of  magnitude)  has  been  subtracted.  The  signals  are  clearly  discernible  until  they  decay  below  10  _5  mV/A,  at  which 
point  they  are  overwhelmed  by  noise. 


where 


Phi, 4)  =  *Hpr(y  -  r ht)] 


(15) 


t=tk 


the  asterisk  denotes  a  convolution  [4],  [42],  and  the  primary  field 
is  that  of  (1)  adapted  to  two  vertically  separated  coils.  An  ap¬ 
proximate  simplified  expression  is  obtained  if  the  time  develop¬ 
ment  of  the  primary  field  is  modeled  by  a  Heaviside  step-off  at 
t  =  0  [43,  p.  118]: 


P(ri,4)  =  hoB (t  =  4)I0]Hpr(r/  -  r ■j,t  =  0 )/J0 

=  Rate  A(4)  R L  Hpr(r/  -  rj)/I0.  (16) 


Our  definition  of  the  time-domain  response  A (tk)  thus  absorbs 
the  medium  permeability  and  Iq,  the  point-to-point  final  value 
of  the  current  that  produces  the  primary  field,  as  discussed  in 
Section  II-A.  A  look  at  Fig.  3  shows  that  a  step-off  response 
like  (16)  represents  a  rather  drastic  simplification  of  what  hap¬ 
pens  in  reality.  Moreover,  we  make  no  attempt  to  model  A  (tk) 
or  express  it  in  terms  of  more  fundamental  quantities.  This  is  a 
difficult  problem  [44]  that  has  been  solved  analytically  only  for 
the  sphere  [40],  [41],  [45]-[47].  We  content  ourselves  with  ex¬ 
tracting  those  principal  elements  and  comparing  them  to  “stan¬ 
dard”  versions  stored  in  a  library  or  verifying  their  case-to-case 
consistency.  The  Appendix  describes  a  fast  and  accurate  numer¬ 
ical  implementation  of  the  model  which  we  use  in  subsequent 
sections. 


B.  Data  Acquisition 

The  MPV  has  been  used  to  take  data  at  several  locations  and 
in  different  experimental  conditions. 

1 )  Vicksburg:  Researchers  from  G&G  Sciences  carried  out  in 
February  and  March  of  2007  a  series  of  test  measurements  [10] 
at  the  Waterways  Experiment  Station  (WES),  a  facility  of  the 
U.S.  Army  Engineer  Research  and  Development  Center  located 
in  Vicksburg,  MS.  They  took  measurements  on  five  targets — a 


Fig.  8.  Measurement  grid  for  the  measurement  carried  out  at  Sky  Hanover. 
Note  that  the  coordinate  axes  are  flipped;  as  explained  in  the  text,  in  this  set  of 
measurements  the  MPV  was  placed  upside-down  and  the  targets  were  moved 
around  on  a  measurement  grid  above.  The  figure  shows  one  typical  setup,  with 
two  40-mm  projectiles  separated  horizontally  by  25  cm. 

BLU-26  bomblet  and  5 7 -mm,  60- mm,  81-mm,  and  105-mm 
munitions — at  different  combinations  of  depth  and  orientation 
using  the  89-point  grid  sketched  in  Fig.  6.  The  test  stand  at 
WES  has  a  machined  fiberglass  platform  over  which  a  sensor, 
mounted  on  a  robotic  arm  with  computer-controlled  motion,  can 
be  positioned  to  an  accuracy  of  1  mm;  the  targets  are  placed  by 
remote  controls,  their  depth  known  to  within  1  cm.  The  back¬ 
ground  was  measured  at  one  faraway  point  before  and  after  each 
data  run.  The  sensor  was  firmly  attached  to  the  robotic  arm  so 
that  its  orientation  was  nominally  constant  throughout.  Fig.  7 
shows  an  example  of  the  raw  (though  background- subtracted) 
data  collected  in  Vicksburg.  In  the  inversions  of  the  next  section 
we  discard  the  first  15  time  channels  of  the  122  that  were  used 
in  this  set  of  measurements.  The  measurements  were  repeated 
27  times  per  data  point  and  thus  lasted  for  2.7  s. 

2)  Sky  Hanover:  Personnel  from  Sky  Research,  Inc.  carried 
out  a  series  of  measurements  at  their  Hanover,  New  Hampshire, 
branch  in  March  of  2009.  This  set  of  measurements  differs  from 
the  others  in  that  the  sensor  was  placed  upside-down  on  an  ad¬ 
justable  support  assembly  and  left  in  place  throughout  the  data 
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Fig.  9.  One  of  the  measurements  performed  in  Vancouver  on  a  cylinder  of  length  12"  and  diameter  3" .  The  MPV  collected  data  on  a  7  x  7  grid;  the  cylinder  lay 
on  a  pivot  of  controllable  angle.  The  experiment  was  performed  at  two  different  target  depths — (a)  and  (b),  respectively.  The  polarizability  elements  are  consistent 
case  to  case.  The  tilt  angles  (measured  values  in  black;  model  predictions  in  shades  of  gray)  are  also  very  close.  Table  I  shows  the  results  of  the  inversion. 


TABLE  I 

Inverted  and  Measured  Depths  and  Directions  of  a  12-in  Cylinder  (a  Small  Component  Along  x  Is  Not  Shown) 


Dip  (°) 

Pmeas 

Pinvr 

^meas  (CHl) 

hmvr  (cm) 

Discr.  (%) 

15 

+0.9  6y  —  0.2  6z 

+0.96*/  —  0.27  z 

41.18 

44.24 

7.43 

30 

+0.87 y  -  0.50 z 

+0.85 y  -  0.52 z 

41.80 

43.26 

3.49 

45 

-0.71*/  +  0.71z 

+0.77*/  -  0.64 z 

50.69 

51.75 

2.09 

60 

—0.50  y  +  0.87  z 

—0.45  y  +  0.89  z 

43.90 

42.66 

2.82 

75 

—0.26*/  +  0.97z 

—0.23  y  +  0.97  z 

44.99 

42.96 

4.51 

15 

+0.9  6*/  —  0.2  6z 

+0.96*/  —  0.2  6z 

25.18 

28.01 

11.24 

30 

+0.87 y  -  0.50 z 

+0.84*/  -  0.54 z 

25.80 

26.78 

3.80 

45 

—0.71*/  +  0.71z 

—0.67  y  +  0.74  z 

27.19 

25.50 

6.22 

60 

—0.50  y  +  0.87  z 

— 0.39*/  +  0.92 z 

27.40 

25.33 

7.55 

75 

—0.25*/  +  0.97z 

—0.19*/  +  0.98z 

28.99 

25.72 

11.28 

collection  process  while  the  targets  were  moved  around,  again 
in  a  snakelike  fashion,  on  a  5  x  5  measurement  grid  with  20  cm 
point-to-point  separation  (see  Fig.  8)  that  rested  on  a  door  atop 
a  stack  of  milk  crates. 

All  measurements  involved  two-target  configurations  and 
employed  30  time  channels  (for  inversion  purposes  we  start 
with  the  second).  Two  stacks,  each  with  27  repetitions,  were 
used,  for  a  measurement  duration  of  5.4  s  per  data  point. 
Some  of  the  targets — the  60-mm,  81 -mm,  and  105 -mm  shells 
from  the  preceding  section — had  already  been  characterized  at 
Vicksburg.  A  40-mm  projectile  was  also  studied  in  this  run,  as 
were  a  ferrous  ellipsoid  and  a  box  of  nails  to  represent  clutter. 

3)  Sky  Vancouver:  Another  set  of  measurements  taken 
by  Sky  Research  staff  took  place  in  their  Vancouver,  British 
Columbia,  office  in  February  and  March  of  2010.  The  objects 
studied  in  this  data  run  were  not  actual  munitions  but  a  diverse 
sampling  of  cylinders  of  sizes  and  material  properties  typical 
of  UXO.  Several  one-  and  two-target  scenarios  were  studied, 
along  with  attitude  configurations  similar  to  the  one  of  Fig.  5 
and  time-on  measurements  to  get  a  complete  picture  of  the 
waveform  (as  presented  in  Fig.  3).  An  assortment  of  grids 
were  used  during  the  measurements,  the  most  common  being  a 
7x7  grid  with  15 -cm  separation,  though  for  the  multi-object 
measurements  the  favored  grid  was  7x6;  unlike  the  previous 
measurements,  in  this  run  the  grids  were  navigated  in  zigzag, 
with  the  instrument  moving  in  the  same  direction  across  every 
transect.  A  typical  measurement  included  27  repetitions  (for 
a  measurement  time  of  2.7  s  per  data  point)  and  contained  35 
time  channels;  in  our  inversions  we  discard  the  first  of  these. 


C.  Results ,  One  Target 

We  first  look  at  some  one-target  inversions  performed  on  the 
Vancouver  MPV  data.  In  one  example,  a  solid  steel  cylinder  of 
length  12"  (30.48  cm)  and  diameter  3  in  (7.62  cm)  was  placed  on 
a  pivot  with  controllable  dip  angle.  (See  the  inset  of  Fig.  9  for  an 
example  run.)  Data  were  collected  on  the  7  x  7  grid  described 
previously  with  the  cylinder  pointing  15°,  30°,  60°,  and  75° 
degrees  below  the  horizontal;  the  experiment  was  performed  at 
two  different  target  depths.  Fig.  9  and  Table  I  display  the  results 
of  the  inversion  for  each  case,  including  the  reconstructed  depth 
and  the  unit  vector  associated  with  the  dominant  eigenvalue;  the 
x  component  is  always  an  order  of  magnitude  smaller  and  is  not 
displayed.  (The  table  also  displays  the  results  of  two  additional 
measurements  carried  out  at  a  45°  dip  angle.)  In  all  cases  the 
agreement  between  measured  and  inverted  locations  and  ori¬ 
entations  is  reasonable  (as  we  can  see  from  the  last  column, 
which  displays  the  percent  discrepancies  in  depth  defined  by 
l^meas  —  ^inv|/^meas  x  100).  The  polarizability  elements  are 
consistent  case  to  case  and  show  the  usual  power-law/exponen- 
tial  decay  [27],  [48];  two  of  the  elements,  the  “transverse”  ones, 
are  very  similar  to  each  other  and  thus  reflect  the  azimuthal  sym¬ 
metry  of  the  object.  Note  that  there  is  an  ambiguity  in  the  orien¬ 
tation  resulting  from  the  cylinder’s  up-down  symmetry. 

We  can  use  the  MPV  data  and  our  dipole  inversion  to  investi¬ 
gate  how  the  material  properties  of  an  object  affect  its  EMI  re¬ 
sponse.  Fig.  10  shows  the  responses  of  two  cylinders  of  identical 
size — length  6  in  and  diameter  3  in — placed  at  similar  “depths.” 
(Two  different  experiments,  at  different  depths  and  dip  angles, 
are  shown.)  One  of  the  cylinders  is  made  of  steel  and  the  other 
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Fig.  10.  Vancouver  measurements  on  (aluminum  and  steel)  cylinders  of  length  6  in  and  diameter  3  in.  The  MPV  collected  data  on  a  7  x  7  grid.  Two  different 
target  depths — (a)  and  (b),  respectively — are  shown  for  each  cylinder.  The  polarizability  elements  of  the  nonferrous  cylinder  are  “inverted”  (in  the  sense  that  the 
longitudinal  element  is  smaller  than  the  transverse  ones).  They  are  also  an  order  of  magnitude  smaller  due  to  the  higher  conductivity  of  aluminum. 


Two  objects: 


—  Pxi  --  Px2 
Pyl  “  ~  Py2 
~  Pzl  Pz2 


Time  (s) 


Time  (s) 


Fig.  11.  Inverted  time-dependent  magnetic  polarizabilities  extracted  from  a 
multi-object  MPV  measurement.  A  105 -mm  UXO  and  a  box  of  nails  are  in¬ 
terrogated  by  the  sensor  on  the  5  x  5  grid  of  Fig.  8.  Two  separate  experiments, 
involving  different  box/UXO  separations  (25  and  40  cm),  are  displayed.  The 
larger  response  is  consistent  with  that  of  the  UXO  in  question;  to  confirm  this 
we  also  plot  the  polarizabilities  inverted  from  two  independent  one-target  mea¬ 
surements  carried  out  three  years  before. 


Fig.  12.  Inverted  time-dependent  magnetic  polarizabilities  extracted  from  a 
multi-object  MPV  measurement.  The  setup  is  identical  to  that  of  Fig.  11,  except 
that  the  UXO  in  this  figure  is  a  60-mm  mortar.  The  independent  one-target  data 
once  again  show  good  agreement  for  the  UXO,  while  the  box  of  nails — shown 
in  the  inset  above — has  a  response  similar  to  that  on  the  other  figure.  Its  quick 
decay  despite  its  significant  early-time  value  can  be  explained  by  EMI  phe¬ 
nomenology. 


is  made  of  aluminum.  Two  features  of  the  polarizability  ele¬ 
ments  of  nonferrous  bodies  stand  out:  1)  The  longitudinal  po¬ 
larizability  element  is  smaller  than  the  two  transverse  ones  for 
a  non-ferrous  body,  the  opposite  of  what  happens  with  a  ferrous 
object  [49];  and  2)  The  polarizability  elements  of  the  aluminum 
cylinder  are  smaller  than  those  of  the  steel  cylinder  by  about 
an  order  of  magnitude.  The  latter  can  be  explained  by  the  fact 
that  the  amplitude  of  the  time  derivative  of  the  magnetic  flux 
through  the  receivers — the  quantity  measured  by  the  MPV — de¬ 
pends  on  the  material  properties  of  the  target  in  the  combination 
[45]  ~a/cr,  where  a  is  a  characteristic  length  and  a  is  the  con¬ 
ductivity.  The  conductivity  of  aluminum  is  a\\  =  36.9  MS/m, 
some  six  times  that  of  steel,  osteei  =  3.9  MS/m  [50].  On  the 
other  hand,  the  higher  conductivity  allows  the  induced  eddy  cur¬ 
rents  to  linger  on  for  a  bit  longer  before  their  eventual  exponen¬ 
tial  decay;  this  is  also  clearly  visible  from  the  figure. 

D.  Results,  Two  Targets 

Figs.  11  and  12  show  the  inverted  polarizabilities  from  some 
of  the  measurements  taken  at  Sky  Hanover  and  described  in 
Section  III-B2.  In  both  cases  a  UXO  is  placed  in  the  field  of 


view  of  the  MPV  alongside  a  full  5-lb  box  of  nails  used  as  a 
stand-in  for  clutter.  In  Fig.  11  the  UXO  is  a  105-mm  howitzer 
shell  and  in  Fig.  12  it  is  a  60-mm  mortar  round;  in  all  cases  the 
munitions  are  horizontal.  Each  figure  shows  two  different  sce¬ 
narios,  with  the  box  of  nails  placed  25  and  40  cm  away  (in  the 
^-direction)  from  the  ordnance  and  at  the  same  depth.  The  fig¬ 
ures  show  the  inverted  time-dependent  polarizabilities. 

In  all  cases  the  “first”  object  (solid  lines)  is  readily  identified 
as  the  corresponding  UXO  from  its  decay  profile.  The  largest 
(“longitudinal”)  polarizability  element  of  the  UXO  is  j3y  in 
every  case,  consistent  with  the  fact  that  the  munition  is  hori¬ 
zontal  and  pointing  in  the  ^rMPV -direction.  (The  corresponding 
eigenvector  is  almost  exactly  ?/MPV  in  all  cases.)  The  trans¬ 
verse  polarizability  elements  of  the  first  object  are  again  very 
similar,  as  expected  from  a  body  with  azimuthal  symmetry. 
To  support  the  foregoing  statements  we  have  included  in  the 
figures  (as  dots)  the  inverted  polarizabilities  of  the  105-mm 
and  60-mm  UXO  in  two  independent  single-object  characteri¬ 
zation  measurements  taken  over  the  89-point  grid  at  Vicksburg 
(Section  III-B1).  The  inverted  locations  are  displayed  and 
compared  to  the  ground  truth  in  Table  II.  (The  systematic  error 
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TABLE  II 

Inverted  Target  Locations  and  Ground  Truth  for  Selected 
Vicksburg  Measurements 


UXO 

Model  (cm) 

x  y  z 

Ground  truth  (cm) 

x  y  z 

105 -mm 

3.94 

3.87 

61.71 

0 

0 

60 

105 -mm 

3.79 

3.40 

42.52 

0 

0 

40 

60-mm 

3.85 

-1.35 

46.71 

0 

0 

45 

60- mm 

3.52 

-1.78 

26.85 

0 

0 

25 

of  almost  4  cm  in  the  rr-coordinate  had  already  been  noted  by 
the  researchers  that  took  the  data  [10].) 

The  “second”  object  (shown  with  dashed  lines)  produces  a 
signal  whose  three  elements  are  clearly  different.  This  is  to  be 
expected,  since  the  box  has  much  less  symmetry  than  the  muni¬ 
tions.  The  polarizability  elements  from  the  second  object  decay 
much  faster  than  those  of  the  first,  even  though  at  early  times 
they  are  of  comparable  magnitude  (and  can  even  be  larger,  as 
in  Fig.  12).  This  behavior  can  be  explained  using  EMI  phe¬ 
nomenology  [48],  [51]:  the  metal  in  the  box  is  packed  in  a  loose 
and  highly  irregular  manner,  occupying  a  fairly  large  but  not 
compact  volume.  (Its  dimensions  are  «  11  x  11  x  14  cm3  and  its 
weight  is  5  lb  or  2.3  kg,  which  makes  it  intermediate  between  the 
two  projectiles.  To  give  a  sense  of  the  packing  fraction  we  note 
that  a  steel  cube  of  those  dimensions  would  weigh  some  30  lb.) 
At  early  times,  when  the  eddy  currents  induced  by  the  sensor 
are  superficial,  the  response  is  sizable;  later  on,  as  the  currents 
diffuse,  they  are  hampered  by  the  relatively  large  hollow  spaces 
between  the  nails  and  dwindle  quickly. 

We  note  that  the  results  are  not  as  consistent  measurement-to- 
measurement  for  the  box  as  they  are  for  the  projectile.  (Note  in 
particular  the  appearance  of  sign  flips.)  The  fact  that  the  box  has 
sharp  corners  makes  it  difficult  to  model  it  using  a  dipole.  Also, 
in  this  experiment  the  sensor  remained  stationary  while  the  tar¬ 
gets  were  moved:  the  separation  between  the  objects  changes 
from  point  to  point,  however  minutely,  and  the  uncertainties  ac¬ 
cumulate;  moreover,  the  nails  within  the  box  also  move,  how¬ 
ever  imperceptibly,  from  point  to  point,  which  also  augments  the 
inversion  error.  When  the  objects  are  placed  at  different  depths, 
the  results  deteriorate  somewhat  for  the  box  but  are  consistent 
for  the  UXO. 

For  a  final  example  we  return  to  the  Vancouver  measure¬ 
ments.  Fig.  13  shows  the  polarizability  elements  extracted  from 
two  different  two-target  experiments.  The  dashed  lines  repre¬ 
sent  “Case  45,”  where  two  copies  of  the  6-in  cylinder  studied 
in  Fig.  10  were  placed  50  cm  apart  at  a  depth  of  52.5  cm. 
The  solid  lines  correspond  to  “Case  67,”  where  one  of  the  6-in 
cylinders  shared  the  field  of  view  of  the  MPV  with  a  smaller 
and  thinner  cylinder  (length  4  in  and  diameter  1.5  in).  For  com¬ 
parison  we  have  included  the  polarizability  elements  of  each 
of  the  cylinders  as  inverted  from  single-target  data  collected 
for  those  objects.  The  measurements  are  again  quite  consis¬ 
tent.  Table  III  compares  the  inverted  object  locations  to  the 
measured  ground  truth.  We  remark  that  there  seems  to  be  a 
small  systematic  error  in  that  the  predicted  depth  is  always 
slightly  larger  than  the  measurements.  Moreover,  some  other 
cases  from  the  Vancouver  run  do  not  lend  themselves  to  easy 
inversion,  especially  for  scenarios  in  which  one  object  is  very 
close  to  another  and  at  a  different  height,  with  the  extreme 


Fig.  13.  Magnetic  polarizabilities  inverted  from  Vancouver  one-target  and 
two-target  measurements.  Case  45  (dashed  lines)  involves  two  copies  of  the 
6  in  x  3  in  cylinder  from  Fig.  10  and  Case  67  (solid  lines)  includes  one  of 
those  cylinders  and  a  smaller  one  of  size  4  in  x  3  in.  For  comparison,  we  have 
added  the  polarizabilities  extracted  from  single-target  data  collected  for  those 
objects.  The  location  predictions  appear  in  Table  III. 


TABLE  III 

Inverted  Locations  and  Ground  Truth  for  Selected  Vancouver 
Two-Target  Measurements 


Model 

(cm) 

Ground  truth  (cm) 

ID 

zi 

Z2 

dz 

dy 

zi 

Z2 

dz 

dy 

45 

56.5 

57.9 

1.4 

50.0 

52.5 

52.5 

0.0 

50.0 

67 

41.8 

34.2 

7.6 

50.0 

38.5 

31.0 

7.5 

50.0 

situation  of  one  target  lying  right  above  the  other  being,  as 
expected,  the  hardest  to  treat  [14]. 

We  have  seen  in  this  section  that  the  MPV  provides  data  of 
sufficient  quality  that  it  allows  the  consistent  identification  of 
deeply  buried  targets.  Using  the  MPV  one  can  also  discern  com¬ 
posite  signals  due  to  two  visually  obscured  objects  in  proximity 
to  one  another. 

IV.  MPV-II 

Some  shortcomings  of  the  original  MPV  were  discovered  at  the 
prototype  stage.  First  and  foremost,  the  sensor  is  somewhat  less 
man-portable  in  practice  than  on  paper.  The  transmitter  and  re¬ 
ceiver  assemblies  have  sizes  and  weights  that  result  in  a  sensor 
head  that  is  simply  too  large  and  heavy.  The  boom  is  attached  to 
the  head  in  a  way  that  places  the  latter  too  far  from  the  operator, 
hampering  the  usability  of  the  sensor  as  much  as  the  fact  that  the 
data  acquisition  system  requires  a  tether.  The  ArcSecond  posi¬ 
tioning  system  described  in  Section  II-C  has  satisfactory  preci¬ 
sion  over  a  short  baseline  but  was  found  to  require  rather  too  much 
attention;  it  failed  to  produce  usable  data  when  one  of  the  three 
AS  I  transmitters  was  obscured,  which  means  that  the  system  may 
lose  much  of  its  usefulness  in  treed  landscapes;  it  also  is  not,  and 
will  likely  never  be,  commercially  available.  Finally,  the  top  re¬ 
ceiver,  R0,  is  always  much  farther  from  any  buried  target  than  the 
other  four  and  of  necessity  detects  a  much  weaker  signal;  the  low 
SNR  may  result  in  data  that  in  some  cases  can  compromise  the 
quality  of  the  rest  of  the  collected  information,  outweighing  the 
gain  in  data  diversity  that  motivated  its  original  inclusion. 

The  MPV-II,  a  new  incarnation  of  the  MPV,  has  been  de¬ 
signed  and  developed  (again  by  G&G  Sciences)  with  these  is¬ 
sues  in  mind;  it  appears  on  Fig.  14.  Its  receiver  cubes  are  smaller 
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Fig.  14.  A  new  incarnation  of  the  MPV,  the  MPV-II,  features  a  smaller  and 
lighter  head,  a  single  transmitter  coil,  and  a  positioning  system  that  uses  the 
transmitter  as  a  beacon.  The  sensor  incorporates  an  enhanced  RTK  GPS  receiver 
and  uses  an  iPad  to  communicate  with  the  data  acquisition  system.  Its  head  is 
about  50  cm  in  diameter. 


(with  8-cm  sides  as  opposed  to  10  cm,  and  with  coil  height  3  cm) 
and  have  been  set  closer  together,  with  a  center-to-center  sepa¬ 
ration  of  18.44  cm  (down  from  39.3  cm).  The  new  system  has 
one  transmitter  coil  instead  of  two,  thus  shedding  much  weight; 
the  21 -turn,  7.74  cm  tall  coil  has  diameter  49.68  cm  (one-third 
smaller  than  the  original)  and  is  centered  with  the  receiver  as¬ 
sembly.  In  all,  the  weight  of  the  new  device  is  about  12  lb,  which 
makes  it  less  than  half  as  heavy  as  the  23-lb  original  MPV.  The 
ergonomic  enhancements  go  beyond  the  reduction  in  size  and 
weight.  The  sensor  head  of  the  new  instrument  is  only  about  a 
foot  away  from  the  operator.  The  setup  is  also  more  balanced, 
for  greater  ease  of  operation.  The  system  uses  an  Apple  iPad  to 
allow  wireless  communication  with  the  data  acquisition  system 
(and  potentially  to  enhance  the  user  interface). 

For  positioning,  the  MPV-II  has  a  real-time-kinematic  (RTK) 
GPS  system  [52]  with  a  horizontal  uncertainty  of  ±2  cm.  The 
device  also  incorporates  an  electronic  compass  that  gives  its 
3D  attitude  at  all  times.  For  enhanced  accuracy,  and  for  situa¬ 
tions  in  which  there  is  no  GPS  access,  a  new  positioning  system 
is  being  implemented  that  uses  the  primary  field  of  the  trans¬ 
mitter  as  a  beacon;  the  field  information  is  collected  by  two  re¬ 
ceivers  separated  by  149.35  cm  and  (at  this  point)  tethered  to 
the  sensor  head.  This  is  similar  to  the  positioning  system  for 
the  GEM-3D+  frequency-domain  sensor  [34]  and  is  depicted 
in  Fig.  14. 

A  final  difference  between  the  MPV  and  the  MPV-II  is  that 
the  zeroth  receiver  has  become  the  fifth;  it  is  now  located  behind 
R1  on  the  same  horizontal  plane  as  the  others.  This  new  config¬ 
uration  of  receivers  should  allow  the  experimental  computation 
of  field  gradients. 

Fig.  15  shows  an  example  dipole  inversion  from  an  MPV-II 
measurement.  The  data  were  collected  by  two  of  the  authors 
(BEB  and  NL)  at  the  Yuma  Proving  Ground,  AZ,  in  October  of 
2010.  A  40-mm  projectile  was  interred  22  cm  below  the  ground, 
vertical  facing  down;  the  MPV-II  was  placed  flat  at  a  set  of  nine 
points  not  on  a  grid  (see  Fig.  14).  The  inverted  depth  and  dip 
angle  are  respectively  22.8  cm  and  268°.  For  comparison  we 
show  a  Sky  Hanover  measurement  (Section  III-B2)  with  two 
copies  of  that  projectile  “buried”  25  cm  away  from  the  sensor 
and  separated  by  40  cm.  The  inverted  locations  are  within  2  cm 
of  the  ground  truth.  The  polarizabilities  are  very  close  in  all  three 


Fig.  15.  Magnetic  polarizabilities  inverted  from  Yuma  Proving  Ground 
MPV-II  data  (markers).  The  target  is  a  40-mm  projectile  buried  22  cm  below 
ground.  The  inverted  location  and  orientation  agree  well  with  the  ground 
truth.  For  comparison  we  include  the  polarizabilities  (solid  and  dashed  lines) 
extracted  from  a  two-target  MPV  measurement  carried  out,  using  the  same 
munition,  at  Sky  Research  in  Hanover. 


cases,  lending  credence  to  both  instruments  and  to  the  inversion 
procedure. 


V.  Conclusion 

In  this  paper  we  have  presented  the  MPV  sensor,  a  time-do- 
main  EMI  instrument  designed  for  the  detection,  identification, 
and  discrimination  of  buried  objects,  in  particular  unexploded 
ordnance.  We  have  seen  that  the  diverse,  high-quality  data  pro¬ 
vided  by  the  MPV  can  be  combined  with  the  tensor  dipole  model 
to  locate  and  characterize  subsurface  targets  accurately.  The 
data  consistently  allow  the  simultaneous  identification  of  mul¬ 
tiple  targets. 

The  newer,  streamlined  MPV-II  is  expected  to  take  over  the 
MPV’s  original  role  in  UXO  remediation,  as  it  produces  data  of 
similar  quality  (and  improving)  and  comes  closer  to  being  truly 
portable  and  usable  in  forbidding  terrains.  The  current  system, 
however,  is  expected  to  keep  being  relevant  in  humanitarian  and 
military-environmental  contexts.  For  example,  it  has  been  pro¬ 
posed  to  use  the  current  MPV  as  part  of  a  robotic  system  that 
could  be  used  to  detect  roadside  bombs  in  theaters  of  conflict 
and,  in  friendlier  contexts,  for  pipe  detection.  In  fact,  the  instru¬ 
ment  has  already  been  used  as  part  of  a  humanitarian  demining 
operation  in  Laos. 


Appendix 

Numerical  Implementation  of  the  Dipole  Model 


As  mentioned  in  Section  II,  at  every  location  r  ?-  and  time  gate 
tk  the  MPV  provides  15  data  points,  corresponding  to  the  flux 
in  the  zvth  direction  through  the  ^th  receiver  R fi.  These  we  call 
tk),  in  analogy  with  (14).  If  measurements  have  been 
taken  at  Npt  different  locations  using  Ntg  time  gates,  we  can 
stack  the  collected  information  in  a  two-dimensional  array  Hsc 
of  size  (TVpt  •  5  •  3)  x  Ntg  with  elements 


Hsc 


(17) 
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where 


Hqv(y\,  t\) 
Hou  (rivpt^i) 

-H±v  (r.v,„  •  /  i  ) 


Hov  (ri,ijvtg) 


H0v  (rNpt,tNte) 
Hu,  (ri ,tNts) 


(18) 


H±v  (rNpt,tNtg)  J 


Knowledge  of  the  location  r'  of  the  responding  dipole  and  of 
the  MPV  geometry  displayed  in  Fig.  2  allows  the  computation 
of  the  (TVpt  •  5  •  3)  x  3  arrays 


with  the  desired  A (tk)  “as  diagonal  as  possible”  [56],  [57].  The 
matrix  V  is  an  approximation  to  Rabc,  which  gives  us  the  orien¬ 
tation  of  the  target. 

As  remarked  right  after  (18),  the  procedure  just  described  as¬ 
sumes  that  we  know  the  location  r'  of  the  dipole.  The  final  step 
in  the  inversion  procedure  is  to  relax  that  condition  and  find  a 
now  unknown  r'  by  optimization.  Given  a  proposed  value  v'p  for 
the  dipole  location  we  can  define  a  “modeled”  signal  by 

S(r£)  =  G(r;)P  (24) 

and  minimize  the  objective  function 


R„  =  [X  Y  Z]  and  Hpr  =  [Hgr  Hpr  Hpr] 

(19) 

along  with  R,  the  row-by-row  2-norm  of  R„,  in  terms  of  which 
we  can  define  the  Green  matrix 


G  = 


(20) 


where 


Gz=  [  ds  -2- 

Jrm  4ttR5 

o,=Ld’  Z? 

G*  =  /a/’sW 


0  [3ZX  3ZY  3Z2  —  R2  ]  0  Hpr 

(21a) 

0[3YX  3Y2  —  R2  3YZ]  0  Hpr 

(21b) 

U  [3X2  —  R2  3XY  3XZ](8)Hpr. 

(21c) 


In  (21),  all  array  operations  (multiplication  and  raising  to  posi¬ 
tive  and  negative  powers)  are  carried  out  element-by-element; 
the  0  operator  is  shorthand  for  elementwise  multiplication  and 
matrix  concatenation  corresponding  to  a  row-by-row  Kronecker 
product  [53].  The  integration  over  the  receiver  areas  is  carried 
out  using  two-dimensional  Gauss-Legendre  quadrature  [54] 
along  a  third  dimension.  As  defined,  G  has  size  (Npt  •  3  •  5)  x  9. 
We  can  then  solve  for  the  polarizability  tensor  using  the  pseu¬ 
doinverse  [55] 

P  =  argmin  1||GM  —  Hsc||2  =  G+Hsc.  (22) 
M  2 

The  resulting  9  x  Ntg  array  P  gives  the  full  polarizability  matrix 
at  all  time  gates.  (It  is  also  possible  to  enforce  the  symmetry  of 
B  and  have  P  have  6  columns  instead  of  9;  the  resulting  Green 
matrix  is  less  transparent,  however,  and  we  did  not  find  a  no¬ 
ticeable  improvement  in  either  speed  or  accuracy.)  We  then  re¬ 
arrange  the  elements  of  P  into  a  three-dimensional  array  of  size 
3  x  3  x  Ntg  and  use  a  joint-approximate-diagonalization  algo¬ 
rithm  that  finds  an  orthogonal  matrix  V  such  that  at  each  time 
gate  tk 


P(4)  =  VA(4)Vt  (23) 


J  = 


\ 


[vec(S(rp  —  Hsc)]Tvec(S(rp  —  Hs 
[vecHsc]TvecHsc 


x  100%, 


(25) 


where  the  vec  operator  turns  a  matrix  into  a  column  vector  by 
stacking  its  columns  in  order  [58].  It  is  not  necessary  to  perform 
any  weighting  or  regularization.  Every  function  call  during  the 
minimization  of  (25)  performs  a  linear  least-squares  fit  of  the 
form  (22),  which  means  that  in  a  sense  each  proposed  dipole  lo¬ 
cation  dictates  at  one  stroke  the  orientation  and  time-dependent 
polarizabilities  that  best  enable  it  to  conform  to  the  measured 
signal. 

Any  unconstrained  optimization  routine  can  in  principle 
be  used  to  extract  the  optimal  r'.  We  have  used  a  subspace 
trust  region  algorithm  based  on  an  interior-reflective  Newton 
method  (implemented  in  Matlab  as  lsqnonlin)  [59]— [61], 
a  Broyden-Fletcher-Goldfarb-Shanno  quasi-Newton  method 
(fminunc)  ([59]  and  references  therein),  the  Nelder-Mead 
downhill  simplex  method  (fminsearch)  [62],  and  differ¬ 
ential  evolution  (DE),  a  genetic-style  minimization  algorithm 
[63].  We  obtain  the  same  results  with  the  different  methods. 
The  procedure  is  also  very  fast — essentially  instantaneous  for 
one-target  inversions. 

It  is  straightforward  to  adapt  the  method  to  two-object  sce¬ 
narios  by  adding  columns  to  the  Green  matrix  and  rows  to  the 
response  array.  Equation  (22)  becomes 


"P1" 

P2 


=  arg^iA  o ll[Gl 

M1  ,M2  Z 


=  [GJ  G2]+Hsc 


G2] 


M1 

M2 


H sc  1 1 2 


(26) 


and  the  minimization  yields  the  two  dipole  locations  iq,  r2. 
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